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Abstract 

We present a determination of \V c b\ from semileptonic B decays that includes resummation 
of supposedly large perturbative corrections, originating from the running of the strong 
coupling. We argue that the low value of the BLM scale found previously for inclusive 
decays is a manifestation of the renormalon divergence of the perturbative series starting 
already in third order. A reliable determination of \V c b\ from inclusive decays is possible if 
one either uses a short- distance b quark mass or eliminates all unphysical mass parameters 
in terms of measured observables, such that all infra-red contributions of order 1 /nib cancel 
explicitly. We find that using the MS running mass significantly reduces the perturbative 
coefficients already in low orders. For a semileptonic branching ratio of 10.9% we obtain 
1^*1(^/1.50 ps) 1 / 2 = 0.041 ± 0.002 from inclusive decays, in good agreement with the value 
extracted from exclusive decays. 
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1 Introduction 



The physics of heavy flavours has experienced a rapid development within the past few years, 
driven by new data that aim to test the Standard Model and to determine its fundamental 
parameters. In particular, semileptonic B decays for the moment provide the best possi- 
bility to determine the CKM matrix element \V c b\. Two competing strategies, which both 
have received considerable attention, are the determination of \V c b\ from the total inclusive 
semileptonic decay rate [IJ and from the exclusive B — > D*W decays at the point of zero recoil 
. In both cases the absence of 1/ rrib corrections allows an accurate theoretical description. 
The decay rates can be calculated within perturbation theory up to terms of order 1/m 2 . 
Moreover, the 1/m 2 corrections are estimated to be rather small (~ 5%). Thus, at present, 
the theoretical accuracy of the determination of \V c b\ is to a large extent limited by a poor 
control over perturbative radiative corrections, which are only known to one-loop accuracy. 
An explicit calculation of the second order correction is a very hard enterprise already for 
b — > u transitions and even more so for b — > c transitions because of the c quark mass, whose 
numerical effect is very important, see || |J. 

A process of major phenomenological interest is the total inclusive B meson decay rate 
with a c quark in the final state, which is calculable in perturbation theory as 



Y{B -> X c ev) = r (a) 



l-C F % (a) + O(al 

71 



where Cp = 4/3 and a = (m c /rrib) 2 - 

The tree-level decay rate, including the phase space factor /i(a), reads 



r fo) 
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la + 8a — a — 12a In a, 
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and the function go(a) is known in analytic form ||. Here and below rrib and m c denote pole 
masses. The perturbative expression in ( |1 . 1| ) should be complemented by non-perturbative 
corrections suppressed by powers of the heavy quark masses ||, and we will take these 
corrections into account in the final analysis. In the major part of the paper, however, we 
restrict ourselves to perturbation theory and estimate higher-order perturbative corrections 
to (0). 

For the realistic value m c /rrib = 0.3, Luke, Savage and Wise |7[ have given an estimate 
for the a 2 correction in ( |1.1|) as 



1.67^-15.1 



7T 



1.3) 



where a s is the MS coupling. For m c /mb — > 0, i.e. for b — > u decays, the estimated size of 
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Figure 1: Generic radiative corrections for heavy-particle decays: (a) leading order, (b) with 
a fermion bubble insertion, and (c) with a chain of fermion bubbles. The broken lines 
represent the lepton pair produced in the decay. 



the second order correction is even more striking 0: 



1-2.41 * v ' -28.7 — . (1.4) 

7T \ 7T / 

The coefficients in front of o? s were in both cases obtained by an explicit calculation of 
the diagrams corresponding to the insertion of a fermion loop into the gluon line in the 
leading-order virtual correction, as in Fig. |l](b), or the splitting of an emitted gluon into a 
light quark-antiquark pair, for the real emission. These contributions are proportional to the 
number of light fermion flavours and the above numerical estimates are obtained by restoring 
the full one-loop QCD /3-function by the substitution Nf — > Nf — 33/2. This replacement 
assumes the hypothesis of BLM [H] that the dominating radiative corrections originate from 
the running of the strong coupling. The result of this procedure is usually expressed as 
a redefinition of the scale of the coupling in the leading-order correction that completely 
absorbs the second order correction. The magnitude of these corrections leads to very low 
BLM scales for semileptonic decays 0: 

fjft™ = 0.07m 6 , ^ c = 0.13m&, (1.5) 

with numerical values of order (350-650) MeV that are hardly acceptable. The authors of 
Ref. 0] interpreted their result as an indication that an accurate determination of \V c b\ from 
inclusive decays requires knowledge of the exact second order correction and even those of 
higher order. In Ref. |§ the large two-loop correction was interpreted as a breakdown of 
perturbation theory which disfavours the inclusive approach to |V^,| in comparison to the 
exclusive one, for which large radiative corrections do not appear in the same approximation 
||. On the other hand, as noted in flip"! , the difference in the size of the a 2 s correction 



for inclusive and exclusive decays largely disappears, when the scale of the leading-order 
correction is chosen equally as ^mcvn^ in both cases. Still, the very fact of low BLM scales 
suggests the investigation of yet higher order radiative corrections. 

It is this question we address in this paper. Our analysis extends the results of Ref. [|7| for 
inclusive decays and repeats that of jllj for exclusive ones in that we resum the effects due 
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to one-loop running of the strong coupling, but to all orders in perturbation theory. Thus, 
our investigation of higher order corrections assumes the dominance of vacuum polarization 
effects also in higher orders and we do not address the question whether knowledge of the ex- 
act two- (and higher) loop corrections as compared to the BLM approximation is important. 
The idea is that if higher order corrections are large at a certain scale, they are presumably 
dominated by running coupling effects and can thus be taken into account exactly, at least 
within the restriction to one-loop running. The remaining corrections are then small and 
therefore can only be accounted for by an exact calculation. Formally, we resum terms of the 
type a s (/?oa s ) n , of which the correction found in 0, is the first term with n = 1. These can 
be traced by a calculation of contributions proportional to given by a chain of fermion 
loops as in Fig. |l|(c). The leading-order BLM scales calculated in |7|, |9| correspond to using 
the QCD coupling at some characteristic virtuality obtained by averaging Ink 2 , where k is 
the gluon momentum, over the leading-order diagram. The resummation that we perform 
in this paper amounts to averaging with the one-loop running coupling a s (k 2 ) itself, rather 
than Ink 2 . We have developed a technique to implement this resummation in Refs. [12, [13 



and refer the reader to these articles for all conceptual and technical issues that we do not 
repeat in the present application to semileptonic B decays. 

We find that the large second order radiative correction to b — > uev transitions calculated 
in is in fact already close to the regime, where the series starts to diverge because of 
factorially growing coefficients. In our approximation (called 'Naive Non-Abelianization' in 
fl2| , p~3|] ) the series in (|1.4j) is continued as 

1 - 2.41^|l + 11.12 + 149.3 + 2319 (^) 3 + 42751 + . . . j, (1.6) 

and with a s {m b ) = 0.21 one gets a non-convergent series of corrections to the decay rate 
already in low orders: 

F(B -> X u ev) = r (0) ( 1 - 2.41 a ^ mb ' [i + 0.75 + 0.67 + 0.70 + 0.87 + 1.27 + . . . 



7T 



r (o){i-2.4i as(mfe) r 
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2.31 ±0.62 



= (0.63 ±0.10) r (0). (1.7) 

In attributing a numerical value to this divergent series we assume that it is asymptotic. Then 
one must truncate it at the minimal term and its value gives an estimate of the intrinsic 
limitation of the perturbative calculation^] which cannot be reduced by computing higher 
orders. 



1 In practice, we adopt a similar procedure based on the Borel integral, and give the principal value of the 
Borel integral as the central value, and the imaginary part (divided by 7r) as an estimate of the uncertainty, 
see Sec. 2 for details. 
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For b — > wez/ transitions the minimal term occurs at third order in a s and its size is 
comparable to the second order correction. Numerically, it gives a 15% uncertainty for the 
total decay rate, not taking into account the uncertainty in the input parameters a s and the 
quark masses. For b —>■ cev transitions, in the same approximation of vacuum polarization 
dominance, we obtain 

F(B -> X c ev) = (0.77 ± 0.05) r (0.3) (1.8) 

with a 7% uncertainty. We also note that the uncertainties associated with the fixed-sign 
factorial divergence of perturbative expansions cannot be reduced by the use of a different 



renormalization scheme, or a change of scale in the coupling ||14|| . Thus, the change of scale 



suggested in Ref. fllQ| , which decreases the second order coefficient, is ineffective already at 
the next order, because the reduction of coefficients is compensated by an increase of a s . 

However, it would be premature to draw a pessimistic conclusion from the apparently 
bad behaviour of perturbative corrections. The large corrections displayed above originate 
from infra-red regions in the integration over loop momenta and produce an uncertainty 
parametrically of order Aqcd/^6- As it turns out, the importance of infra-red regions 
is solely due to the choice of an input parameter - the pole mass as renormalized mass 
parameter -, which is incompatible with the short-distance properties of the decay process. 
The series that relates the pole to the bare mass contains large finite renormalizations of 
infra-red origin. If these are made explicit - for example using the MS renormalized mass 
- they cancel with the large corrections of infra-red origin present in the perturbative series 



for the decay width [15, 16 



The preference of the MS (or another 'short- distance') mass might seem surprising and 
even counter-intuitive. After all, it is the pole mass that governs the (partonic) decay 
kinematics and it is the visualization of an almost on-shell (up to effects of order Aqcd) b 
quark inside the meson that motivated the approximation of the meson decay by a free quark 
decay in the first place. But this picture also implies the existence of a static field (in the rest 
frame of the quark) around the heavy quark, which behaves as 1/r at short distances. Thus 
a contribution of order Aqcd to the self-energy of the quark is stored at large distances, r ~ 
1/Aqcd, of the order of the radius of the heavy-light meson. However, due to the Kinoshita- 
Lee-Nauenberg cancellations, in an inclusive decay of a heavy quark the decay vertex is 
localized to within a distance l/m& and the energy stored in the field at large distances 
r > l/n (where > fi > Aqcd is a factorization scale) cannot participate in the hard 
process, but is absorbed into a rearrangement of the colour field of the hadronizing spectator 
quark. This explains, loosely spoken, why a short- distance mass is more appropriate in the 
description of inclusive decays as a hard process. This reasoning assumes that the quark 
produced by the weak current is fast and does not apply to a massive quark produced with 
zero recoil (cf. Sec. 3.2). 

Since ultimately any quark mass parameter is unphysical, the most transparent way 
to exhibit the infra-red cancellations would be to eliminate any mass parameter in terms 
of a suitable physical quantity, provided it is determined by short distances and does not 
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import l/m& corrections (which rules out meson masses). However, as with the strong 
coupling constant, it is convenient to use a mass parameter for book-keeping purposes and 
as a numerical input parameter. It is only at this point that the MS mass (or any other 
mass defined at short distances) is favoured over the pole mass, because it does not import 
infra-red effects. 

In fact, the divergence of the series of radiative corrections to the decay rates is only 
one aspect of the problem with using the pole mass. Another aspect is that it cannot be 



accurately extracted from measurable quantities through perturbative expansions [[L5|, [L7| . 

To reinforce this point we imagine that we used pole masses as numerical input param- 
eters, determined from another measurement. Then one would always find that the size of 
perturbative corrections does not allow a determination of the masses to an accuracy better 
than ±100 MeV (we quote the estimate from |13]|). Treating the uncertainty in the input pa- 



rameters as uncorrelated with the uncertainty in the theoretical prediction for the radiative 
corrections to the B decay width, we obtain a ±21% and ±10% uncertainty for B — > X u ev 
and B — > X c ev decays, respectively. This translates into an irreducible theoretical uncer- 
tainty of ±10% in \V u b\ and ±5% in \V c b\. The 10% uncertainty in T(B — > X c ev) is to be 
compared with the present experimental uncertainty of 12% of the branching ratio. 



The calculation of decay rates with running masses has already been considered in [18| at 
the level of one-loop perturbative corrections. After resummation, we find that apart from 
the cancellation of infra-red contributions the perturbative coefficients are strongly reduced 
already in low orders. In particular, if the MS running mass is used in the tree-level decay 
rate, the series of radiative corrections for b — > uev decays becomes 



1 ±4.25-^jl± 8.99 \-^J ±35.15 + 241.1 + 1547 ±...j. (1.9) 

Although the leading-order correction has increased (and changed sign), the higher-order 
coefficients are significantly reduced, so that with the same value of a s as above we get the 
well convergent series 



7T 



1 ± 0.604 ± 0.159 ± 0.073 ± 0.032 ± 0.022 ± . . . 



r (0) { 1 + 4.25 tts(m ^ \l.92 ± O.Oll 1 (1.10) 



7T 



with an uncertainty that is negligible compared to the uncertainty inherent in the restriction 
to vacuum polarization corrections. 

It should be noted that anticipating the eventual cancellation of infra-red regions we can 
define the numerical value of the sum of radiative corrections even when it diverges, with 
some (ad hoc) prescription. Provided we use the same prescription to define the pole mass, 
and provided we know that the cancellation occurs, we can then simply delete nearly all 
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large uncertainties. Thus, the calculation in the on-shell scheme (using pole masses) can be 
saved at the cost of introducing consistent non-perturbative prescriptions to define the pole 
mass and to sum the series of radiative corrections. We shall consider this (conceptually less 
appealing) possibility as well, and shall see that after resummation of (3q corrections 
the decay rates calculated in the on-shell and MS schemes are close to each other, provided 
the same short- distance mass is used as an input parameter. This supports a posteriori the 
assumption that the dominant higher order corrections are taken into account by vacuum 
polarization effects. 

We conclude that the large corrections found in |7j do not endanger the accuracy of the 
theoretical treatment of inclusive decays. We do find, however, large corrections beyond 
second order in a s , in particular in the on-shell scheme, defined in the sense of the previous 
paragraph. We suggest that these corrections are more important than the a 2 s corrections 
left out by the restriction to the effects of running coupling. Although the accuracy of this 
restriction is not known and is certainly the main deficiency of our analysis, we believe that 
resummation of one-loop running effects provides a fair estimate of higher order perturbative 
corrections and the corresponding 'M-factors' (defined below) should be incorporated into 
any phenomenological analysis. 

As already emphasized above, an appropriate treatment of quark masses as input pa- 
rameters is of equal importance as the size of radiative corrections to the width itself. Most 
previous determinations of | V^^ | from inclusive decays have used the OS scheme and pole 
masses as numerical input. Instead, we choose the MS masses as numerical input parame- 
ters. Thus, when we use the OS scheme in the above sense, the pole masses are calculated 
from the MS masses. This procedure is not without its own difficulties, since we must rely 
on direct determinations of MS masses, as from QCD sum rules, which have been obtained 
without the resummation which we implement for the decay width. 

We use our results for a new determination of \V c b\ from both inclusive and exclusive 
decays. We find a good agreement between the determination in the MS and OS scheme, 
and obtain as our final value 

|Kb|(Wl.50ps) 1/2 = 0.041 ± 0.002, (1.11) 

where the combined error comes from several sources that will be detailed below. 

The presentation is organized as follows. In Sec. 2 we review the necessary formulas 
for resummation. Section 3 contains our main results for the BLM-improved perturbative 
series in B meson semileptonic decays. These two sections are more technical and those 
readers interested only in results may continue with Sec. 4 directly. The updated analysis 
of \Vcb\ obtained with the resummed formula along with its major uncertainties is given in 
Sec. 4, while Sec. 5 is reserved for a summary and conclusions. Some technical discussion 
and especially long formulas are given in the appendices. As a new analytic result, we derive 
an expression for the total b —>■ u semileptonic decay width with a non-zero gluon mass, 



which provides the input necessary for resummation of running coupling effects [12, 13] 
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2 General formulas 

We now formulate the resummation in precise terms. We are interested in the effect of 
radiative corrections to the total semileptonic width, which we define as 



t(b - x c eu) = r (i - r, ^4„ 



7T 



Ed 

n=l 



(2.1) 



The functions d n (a) depend on the ratio of the charm and bottom quark masses, and are 
polynomials in the number of light flavours Nf-. 



d n (a) = d n0 + d nl N f + . . . + d nn N r f 



(2.2) 



The coefficient d nn (a) comes from the insertion of n fermion loops in the gluon lines in the 
leading-order correction; this is the quantity we calculate explicitly. Substituting Nf — > 
Nf — 33/2 in the highest power of Nf, we rewrite Q2.2| ) as 

d n (a) = 8 n (a) + (-p ) n d n (a) , (2.3) 

where the (uncalculated) 5 n is by construction at most of order NJ~ . We use the definition 
of the first coefficient /3q of the QCD /3-function including the factor — 1/(4-71"), 



A = _-L(ii_| W/ 



(2.4) 



with Nf = 4 for the case of interest^. The standard BLM prescription || uses d\ to fix 
the scale in the coupling in the leading-order correction. In the generalization developed in 



Y\ p4| [13|1 all terms 5 n in (|2.3| ) are neglected, but all d n are kept as an estimate of radiative 
corrections for arbitrary n. Then we get 

r(B^^) = r (i-c i , a ' (m » ) 



7T 



g (a) l + J2(-Po) n d n (a)a^(m b ) 



n=l 



(2.5) 



Note that because of the factor l/(47r) in j3 the expansion parameter is effectively a s /(47r). 
To quantify the effect of partial summation of N orders, we introduce the 'M-factors' 



A? 



M^p c [a, -(3 a s {m b )} = 1 + £ {-(3 ) n d n {a)a n s {m b ) 



71=1 



M^ c [a,-f3 a s (m b )} = MfcZJa, -f3 a s (m b )] 



(2.6) 



that measure the modification of the leading-order radiative correction by integrating with 
the running coupling at the vertex. The limit iV — * oo in Eq. (|2.6| ) does not exist in a rigorous 



2 We neglect the charm quark mass in quark loops. Its effect is small (Tj 
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sense, reflecting the factorial divergence of the coefficients d n in high orders. Assuming that 
the perturbative series is asymptotic, one is led to the conclusion that the uncertainty in the 
summation is in fact power-suppressed in mf,, and can be estimated numerically. Thus, in 
the following, numerical values of M ro will always be given with an uncertainty, reflecting 
this problem. This uncertainty cannot be eliminated without a rigorous factorization of 
the corresponding infra-red contributions into the matrix elements of higher dimensional 
operators. 

An equivalent way to present the results is to absorb the M-factors into a redefinition of 
the scale in the lowest-order correction 3 



N 



a s (m b )M^ c [a, -(3 a s (m b )] 



H'oo 



b^c 



(2.7) 



The scale fi\~^ c is just the leading-order BLM scale studied in |7j and the \x n with n > 1 
correspond to a more accurate treatment of the distribution in the gluon virtuality, reflected 
by the size of higher-order corrections with up to n fermion loops. The uncertainty in the 
summation of the series is translated to the uncertainty in the ultimate BLM scale /ioo [12|. 

The calculation of the coefficents d n (a) requires the evaluation of diagrams such as those 
in Fig. [1], with the insertion of n fermion loops in the gluon lines. This problem is solved in 
a most economical way by applying a dispersion technique, and reduces to the calculation 
of the leading-order diagrams with finite gluon mass A. Denote by 

-T (a)C F a s /7r g (a)d (a, A 2 ) 

the sum of the diagrams in Fig. || calculated with a finite gluon mass A, so that do(a, A 2 = 
0) = 1. For the contribution from the one fermion loop insertion, Smith and Voloshin [B5| 
have derived a useful representation (in the V scheme of H) 



^) = -/ o °°^U(a,A 2 )- 



A 2 + ml 



(2.8) 



which has been used in the analysis of Ref. [[?[]. The calculation of diagrams with multiple 
fermion loop insertions involves precisely the same function d (a, A 2 ), and can thus be done 
at little additional calculational expense. In particular, the fixed-order coefficients d n (a) are 
RBI 



obtained as 



d n (a) 



d n 
du n 



B[D](a,u) lu=0 , 



B[D](a,u) 



sin(7ra) r°° d\ 2 ( A 2 



7T 



o A 2 \/i 2 



dnia. A ) — 1 



For finite N one must expand (2.7) in a s and truncate the expansion at the desired order. 
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II 





III 



IV 



Figure 2: Leading-order radiative corrections to the transition operator, whose imaginary 
part gives the inclusive semileptonic decay width. Double line: b quark, solid line: c quark, 
dashed lines: leptons. 



ITU JO \fl z J 



(2.9) 



where d' (a, A 2 ) = (d/d\ 2 )d (a, A 2 ) and C is a scheme-dependent finite renormalization con- 
stant. In the MS scheme one has C = —5/3, in the V scheme C — 0. It is easy to check that 
for n — 1 the above expression reproduces Eq. fl2.8|) . 

A closed expression can be derived for the sum of all diagrams with an arbitrary number 
of fermion bubble insertions [12, 13]: 

(-(3 a s )M oo [a,-(3 a s ] = / d\ 2 $(A 2 ) d' Q (a, A 2 ) + [d (a,X 2 L ) - 1], 

Jo 

where a s = a s (/x), 



$(A 2 ) = arctan 

7T 



l-(3 a s \n(\ 2 /n 2 e c ) 



A 2 ) 



and 



A 2 - - 2 - 



(2.10) 

(2.11) 
(2.12) 



^ L - -fi exp[l/((3 a s ) - C] 

is the position of the Landau pole in the strong coupling. Note that the term with the 9 
function exactly cancels the jump of the arctan at A 2 = — X 2 L , so that <3>(A 2 ) is a continuous 
function of A 2 . 

In this paper we cannot give a detailed discussion of the assumptions underlying the 
derivation of Eqs. (|2.9| ) and fl2.10| ), and refer the reader to the corresponding sections in 
]. Still, two short comments are appropriate. 
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First, note that the product a s (/x)Moo[a, — /3oa s (//)] is explicitly scale-invariant, provided 
the running of the coupling is implemented to leading logarithmic accuracy: a s {fi\) = 
cn s (l^2) /(I — /9o«s(A i 2) hi(/i 2 / /i 2 ))- This result is also scheme-invariant, and in particular inde- 
pendent of the renormalization constant C, provided the couplings are consistently related 
in the same BLM approximation, that is by keeping only the terms with highest power in 
Nf. This is in contrast to the finite order summation coefficients Mjy, which are scheme- 
and scale-dependent. In the following we assume the MS scheme for the coupling a s , and 
the normalization point /i = mj,. 

Secondly, notice that the second term in ( |2.10| ) involves the radiative correction to the 
decay rate with a finite gluon mass, analytically continued to the Landau pole A 2 < 0. The 
renormalon divergence of the perturbation theory is reflected |I6| by non-analytic terms in 
the expansion of do(a, A 2 ) at small A 2 and leads to an imaginary part in this continuation. The 
size of the imaginary part (divided by tt), SM^ = l/(7r|/3o|a a ) lmdo(a, A|), yields an estimate 
of the ultimate accuracy of perturbation theory, beyond which it has to be complemented 
by non-perturbative corrections. The real part of (|2.10|) coincides with the sum of the 
perturbative series defined by the principal value of the Borel integral | 13[] , and the imaginary 
part of do(a, X 2 L ) coincides with the imaginary part of the Borel integral. 

The calculation of the diagrams in Fig. |2] with a finite gluon mass is straightforward, 
albeit tedious, and has been undertaken in 0. Since no formulas were given there, we had 
to redo this calculation. For b — > uev decays, that is for a massless quark in the final state, 
we have succeeded in obtaining an analytic expression for the decay rate. For b — > cev 
decays, we leave the answer in form of at most two-dimensional integrals and evaluate them 
numerically. The corresponding formulas are collected in Appendix A. 



3 Resummation of BLM-type radiative corrections 

Our results are summarized in Table [l|. For several representative values of ^Ja = m c j we 
give calculated values of the fixed-order coefficients d n , the partial sums M n and the BLM 
scales /i n . For n = 1, the /3 a; 2 correction, our results coincide with the ones obtained in J7|. 
We now discuss the numbers in Table [TJ in detail. 

3.1 Hierarchy of BLM scales 

The BLM scale /ioo can be larger than the leading-order scale fi\. This may come unexpect- 
edly. 

The (leading-order) BLM prescription uses the average of Ink 2 , the average virtuality of 
the gluon, as the scale in the coupling. In high orders, this substitution generates a series of 
radiative corrections with a geometric growth of coefficients d n ~ {di) n to be compared with 
the factorial growth of the exact coefficients. Resummation to all orders corrects for this 
discrepancy by adjusting fi^ so that the expansion of a s (/Xoo) gives the correct d n for all n. 
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y/a 





0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.99 


9o 


1.81 


1.63 


1.42 


1.25 


1.12 


1.01 


0.92 


0.85 


0.75 


d 


1 


1 


1 


1 


1 


1 


1 


1 


1 


d x 


5.34 


5.00 


4.55 


4.09 


3.59 


3.10 


2.62 


2.16 


1.18 


d 2 


34.4 


30.9 


27.3 


23.9 


20.5 


17.4 


14.3 


11.4 


5.24 


M 


1 


1 


1 


1 


1 


1 


1 


1 


1 


Mi 


1.75 


1.70 


1.64 


1.57 


1.50 


1.43 


1.37 


1.30 


1.17 


M 2 


2.42 


2.31 


2.17 


2.04 


1.90 


1.77 


1.65 


1.53 


1.27 


Moo 


2.31 


2.35 


2.24 


2.10 


1.96 


1.82 


1.69 


1.57 


1.38 




±0.62 


±0.52 


±0.44 


±0.37 


±0.32 


±0.26 


±0.21 


±0.15 


±0.03 




0.07 


0.08 


0.10 


0.13 


0.17 


0.21 


0.27 


0.34 


0.55 




0.13 


0.13 


0.14 


0.15 


0.17 


0.20 


0.23 


0.27 


0.37 



Table 1: Resummation of /3qQ;™ +1 corrections for semileptonic B decay widths, see text. The 
values of M n and fi n are given for —ffi a s {mb) = 0.14. 



Since for this reason for large n the true d n will always outgrow {di) n and since the series is 
with fixed sign, one might suspect that the usual BLM scale-setting rather underestimates 
higher-order corrections. 

For small c quark masses the effect is opposite, and the very small leading BLM scale fii 
is simply an artefact of truncating the perturbative expansion at low order. The scales \i\ 
and /^oo are given by 



/ii = nib exp[— d\/2] = nib exp 



2(3 a ; 



\M 1 



Hoc = m b exp 



2p a s 



1 - 



Mr. 



(3.1) 



Although SM = Mqo — Mi is positive in all cases we consider, the expression in parentheses 
in the first equation in ( |3.1|) is numerically larger than the similar expression in the second 
equation, provided 5M < (Mi — l) 2 /(2 — Mi). This is satisfied for small mass ratios, see 
Table [l], and results in /jl^ > \i\ (recall that (3q is negative with our definition). 

Note that the scale /ioo is bounded from below: as long as M m is positive, /Xqq is larger 
than nib exp[l/(/3oa: s (m&))], the position of the infra-red Landau pole in the running coupling. 
There is no such restriction for fi±, which can take values below the pole. Again, this is an 
artefact of the truncation at fixed order. 
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3.2 Suppression of infra-red contributions for finite c quark mass 



With a massive quark in the final state the radiative corrections apparently are reduced: 
The M-factors M n are smaller, and the BLM scales fi n are larger. This has been observed in 
|7j for the BLM scale in leading order, and it continues to all orders, although the difference 
between m c /mb = (relevant to b — ► uev transitions) and m c /mb = 0.3 is less pronounced 
after resummation. The ambiguities related to the summation of a divergent series are also 
reduced and almost vanish in the limit of zero recoil m c —> m b . One can understand this by 
continuing the argument given in the introduction. As explained, these ambiguities arise, 
because the use of the pole mass parameter implies a static picture, while the energy stored 
in the field at large distances cannot be converted into hard radiation in the weak decay 
process. This assumed that the produced quark is fast in the rest frame of the initial b 
quark. When m c — > m^, the c quark is slow. Then, since the long-range part of the field of 
the b quark is universal, it can be smoothly transferred to the c quark and is simply irrelevant 
for the description of the decay. Therefore these long-distance contributions cannot be seen 
in the form of ambiguities in the zero- velocity limit. 

To see this more explicitly, we recall that contributions of small momenta to decay rates 
can be traced by non-analytic terms in the expansion at small values of the gluon mass [|16 |. 
For the leading-order radiative correction to the B decay width this expansion takes the 
formf] 



do(a, A ) = 1 + hi(a] 

with 

hi(a) = 



\ 



A 2 



nn, 



A 2 

+ h 2 (a) — k + h 3 i(a) m(A 2 /m?) + h 32 (a 
mi l 



A 



2 \ 3/2 



+ 0(A 4 lnA 2> 



mt 



(3.2) 



7T 



2/i(a)#o(a 
-48a 3/2 In a - 12a 2 In a 



5 - 16a 1/2 - 24a - 24a 3/2 + 24a 2 + 48a 5/2 - 8a 3 - 8a 7/2 + 3a 4 

(3.3) 



Note that the tree- level phase space factor fi(a) and the leading-order radiative correction 
go(a) are extracted. The function h\(a) is plotted as a function of the mass ratio \fa = m c /rrib 
in Fig. H For the realistic value m c /mb = 0.3 it is reduced by approximately a factor 2 
compared to the massless case. 

These long-range contributions to the static field are responsible for the major part of 
the ambiguity in the sum of the perturbative series. Within our approach, this ambiguity 
is related to the imaginary part of do(a, A 2 ), continued analytically to the position of the 

4 Note that terms 0(A 2 In A 2 ) are absent ||l6| . This is explained by the absence of a renormalon ambiguity 
in the kinetic energy of the heavy quark inside the B meson, at least in the approximation considered here. 
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hi 




0.4 0.6 

m c /m b 



Figure 3: Coulombic contributions to the B decay rate, Eq. ( |3.3| ), expressed in terms of pole 
masses, as a function of the mass ratio m c /m b . 



Landau pole (|2.10|) , (|2.12|) , and equals 



A/ 



-Ajc^Moo = hx{a)— + 0(\t/mi). 
m b 



(3.4) 



The decrease of the value for at y/a = 0.3 {b — > cev decays) compared to a = {b — > uev 
decays) roughly equals the decrease of the uncertainty. 

In fact, these infra-red contributions are spurious, and can be removed by re-expressing 
the decay widths in terms of the short-distance (say, MS) b and c quark masses instead of 



the pole masses [H], [Rj|. To trace this cancellation we write e.g. the b quark pole mass m b 
as related to the running MS mass nib by the perturbative expansion 



m„ - 777,,(777,,)0 : C F ^^r (A 2 ) + ...), 



7T 



(3.5) 



where we keep a finite gluon mass A as for the decay width. The expansion of r (A 2 ) at small 
gluon masses reads |T5|, |TE| |: 

ttA 



r (A^ 



1 



2rrih 



(3.6) 



Using Eqs. (|3.5| ), (|3.6| ) and similar expressions for the c quark, we find that, when calculated 
with a small gluon mass, the tree level decay rate is modified to 



r (a,A) = r (a){l-CV^ — 

7r m b 



5/2 + (^-a) — lnA(a) 
da 



+ ... 



(3.7) 



It is easy to see that the correction linear in A exactly cancels with a similar term in the 
radiative correction to the decay width, 



r = r (a,A) 



Qi 

1 - C F — g (a)d (a, A 2 ) + 

7T 
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n 


r n 


M b 

n 


d n (0) 


M b->u 

n 


d n (0) 


M b u 

n 





1 


1 


1 


1 


1 


1 


1 


4.6861511 


1.656 


5.3381702 


1.747 


4.3163 


1.604 


2 


17.622650 


2.001 


34.409913 


2.422 


8.0992 


1.763 


3 


109.85885 


2.303 


256.48081 


3.126 


26.680 


1.836 


4 


873.92393 


2.638 


2269.4131 


3.997 


82.262 


1.868 


5 


8839.6860 


3.114 


23679.005 


5.271 


421.33 


1.890 


6 


105814.28 


3.911 


289417.40 


7.450 


1656.1 


1.903 


7 


1484968.4 


5.476 


4081180.2 


11.75 


12135 


1.916 


8 


23740736. 


8.978 


65496131. 


21.42 


52862 


1.924 


oo 




2.066 ±0.231 




2.314 ±0.615 




1.925 ±0.012 



Table 2: Effect of the elimination of the b quark pole mass on the radiative corrections to 
b — ► uev decays, see text. The given values of M n correspond to — (3q^ a s {m b ) = 0.14. 



r (o) 



c, 



a s 7rA 
tt nib 



|5/2 + (Va- 



aj — ln/i(a) ± -g (a)hi(a) \ + ... 



da 



71 



(3. 



with the function h\{a) defined as above. The terms in curly brackets add to zero, so that 
the total decay rate is free from infra-red contributions to this accuracy [15, 16]. It is only 



the use of a pole mass as an input parameter which introduces infra-red l/rrib effects in the 
tree-level decay rate and in the radiative corrections, which cancel in the product. For a 
massless quark in the final state several terms in the small-A expansion can easily be obtained 
analytically, with the result 



r = r (o) li+c F 



0:.s 



7T 



65 




+ 47rln(A 2 /m2))^ 
' mi 



+ . 



. (3.9) 



The infra-red contributions now start at order 0(X 3 /ml In A 2 ) and are numerically negligible^. 
If so, it is natural to formulate the perturbative calculation in terms of a mass parameter 
defined at short distances, so that large infra-red contributions do not appear. We address 
this task now. 



20 



The remaining (small) uncertainty is related to contributions of dimension 6 operators to the decay rate, 
which produce non-perturbative c orre ctions of order 1/ro 3 . They can be relevant for D decays, see 
The terms proportional to A 3 in (3J3) produce an uncertainty of order 10% in the decay rate D 
which can be taken as an indication of the minimal size of 1/rn 3 corrections. 



3- 

Xev, 
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3.3 Elimination of the pole mass: b —> ueu decays. 

The b quark pole mass is related to the MS mass by the perturbative series 

a s {rn b ) r 



m b = m b (m b ){ 1 + C F - 



7T 



1 + ]T f n a n 8 ( 



m) 



n=l 



As above, we approximate 



(3.10) 



(3.11) 



where r n corresponds to contributions of n fermion loops to the leading-order diagram for 
the fermion self-energy and can be calculated using a representation similar to (|2.9| ), ( |2.10|) 
in terms of the leading-order diagram with a finite gluon mass ro(A 2 ), see [12 , |13|1 for details. 
The partial sums for the perturbative series truncated at order N are defined as 



A' 



M b N [-p a s (m b )} = l + ^2(-Po) n r n a:(m b ) 



n=l 



M£,[-/3 a.(m 6 )] = M 



-(3 a s (m b )} . 



(3.12) 



The coefficients r n were calculated in Ref. []T^| and are given together with the partial sums 
in the second and third columns in Table |2|. The perturbative series defining the pole 
mass is divergent [15|, [17], which is reflected by the uncertainty in the factors M^. The 
crucial point is that these uncertainties in defining resummed pole masses are correlated with 
uncertainties in the resummed radiative correction to the decay rate, and cancel against each 
other, when the pole mass is defined by its relation to the short-distance mass as in Q3.10 ) 
or eliminated in favour of the MS mass [pl|, [16| . 

In what follows we shall consider both possibilities. The first one, which we refer to as 
calculation in the on-shell scheme (OS), is to define the resummed inclusive decay rate as 



T os (B -> X u ev) = V, 



OS 



a 



C F ^g (0)M' 



71 



b—tu 
oo 



(3.13) 



where it is understood that the b quark pole mass appearing in the tree-level decay rate is 
substituted by 

(3.14) 



m b = m b [m b ) 



a. 



1 + C F —M: 



7T 



with the factors M^ u and given in Table |2| and the uncertainties deleted. 

The second possibility is to use the MS scheme from the very beginning. Using ( |3.10| ) we 
can write the decay rate as 



T(B -> X u ev) 



7T 



JV 

i + E( 

n=l 



-f3 ) n d n (0)a:(m b ) 



(3.15) 
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where 



3o(O) = 5-0„(O), d n (0) 



go(OK(Q) - 5r r 



(3.16) 



and the tree-level decay rate T Q is expressed in terms of m(,(rn&). Note that the leading-order 
radiative correction changes sign and becomes somewhat larger. 

The coefficients d n and partial sums M N u defined in an obvious way in analogy to 
( p.6| ) are given in Table in comparison to d n and M^ u , respectively. It is seen that the 
coefficients are drastically reduced, and the series has become well convergent. The remaining 
infra-red effects, relevant to the divergence of the perturbative series, are suppressed by three 



powers of the b quark mass |L6| and have become tiny. Most importantly, this improvement 
is effective already at n = 2. We conclude that perturbative coefficients in the MS scheme are 
likely to be much smaller than in the OS scheme without restriction to vacuum polarization 
corrections, too. 

In the framework of a purely perturbative calculation the use of the OS scheme can only 
be justified up to the order where perturbative series diverge. In all-order resummations 
such as the one considered in this paper, one must make sure that the prescription defining 
the pole mass in terms of a short-distance mass or any physical quantity is consistent with 
the prescription to sum the perturbative series of radiative corrections to the decay width. 
Even in this case, the OS scheme is somewhat unnatural since it involves large cancellations 
between radiative corrections to decay rates and to the pole masses already in low orders. 



The resummed decay rate in the MS scheme is readily obtained by inserting ( |3.14|) into 
( |3.13| ) and expanding up to 0(a s ): 



O'.s 



r n + c F — 



71 



5M» 



<?o(0)M ( 



(3.17) 



The difference between ( |3.13| ) and ( |3.17| ) is an effect of order (CVa s /7r) 2 , which is beyond 
our accuracy. It is a pure scheme dependence, resulting from our incomplete perturbative 
calculation. Numerically, the difference amounts to about 6%, which is significantly smaller 
than the ±15% uncertainty for the radiative corrections noted in the introduction. 



3.4 Elimination of the pole mass: b — ► ceu decays. 

Expressing the h — > cev decay rate in terms of the running masses is slightly more cumber- 
some. As above, we start with the resummed decay rate in the OS scheme 



T os (B -> X c ev) = T 



OS 



C F ^g (a)M ( 



7T 



oo 



(3.18) 



where now both the c and b quark masses have to be expressed in terms of the running 
masses. Here we want to be somewhat more general, and introduce running masses at 
arbitrary scale ji as 

m c>b = mM [l + Cf^^M^)] . (3.19) 



7T 
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The factors M^(/i) and M oc (/i) can most easily be calculated by observing that the pole 
mass on the left-hand side of (|3.19|) is scheme-invariant, and using the renormalization- group 
expression for the running mass 



m(n) = m(m) exp f --j- f ^7 m (a) ) 

\ PO Ja a {m) a A J 



(3.20) 



where to our accuracy we have to expand the exponential to first order, but keep the mass 
anomalous dimension to all orders in /3q«" +1 , which is known from [22], [13j (see Eq. (4.19) 
in 



my- 



7T 



5 35 
l--((3 a s )- — ((3 a s ) 2 + ... 
o 3d 



Thus, we obtain, e.g. for the c quark 



7T n /Jo Ja a (m c ) CT 



da 



(3.21) 



(3.22) 



where is a factor relating the c quark pole mass to the running mass m c (m c ), defined 
as in ( pn4|) . 

Collecting everything, we find 

CK s (/i) 



r(£^X c eP) = r (a,//) 1 + C F - 



7T 



5M oo 0i)- Sb (o)Mi 



oo 



+ 2(M C 00 (/i)-M 00 ( / ,))a^ln/ 1 (a) 



(3.23) 



where a = (m c (/z)/m b (/z)) 2 . It is worth while to note that the factor that appears in front 
of the logarithmic derivative of the phase space function f\ is scale-independent: 

7T 



a s (m b ) da 
PoL'F Ja a {m c ) a 

The expression in ( |3.23|) is to be compared with the leading-order decay rate [[18 



(3.24) 



7T 



g (a,n) 



9oW 



15, m? 3 
In — 



a In a— In /i(a) . 



(3.25) 



4 ~" /i 2 2 (ia 

To quantify the effect of resummation, we introduce the corresponding M-factor by rewrit- 
ing the resummed result in ( [3.23j ) as 



T(B - X c ev) = r (a,/i) {l + C F ^-g (d, /i)M^ c (/i) 



7T 



(3.26) 



We give the corresponding values in Table [7] for /i = m b . Tracing the scale dependence of 
Moo(/i) is misleading in this case because g (a, fi) is strongly scale-dependent. 
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3.5 The SV limit and exclusive B decays 



The inclusive decay rate in the Shifman-Voloshin (SV) limit m&, m c ^> m b — m c ^> Aqcd 
p3| is dominated by two exclusive decay channels, B — > Dev and B — > D*ev. Since the final 
state meson is produced almost at rest (in leading order in the heavy quark expansion), the 
techniques of the Heavy Quark Effective Theory are applicable, yielding the decay rate 



T(B -> Deu,D*eu) 



G 2 F \V cb \ 2 (m b 



60tt 3 



(3.27) 



where rjv and t]a are the short- distance matching coefficients of the QCD heavy-heavy cur- 
rents to the corresponding currents in the effective theory (at zero recoil): 



cj/J> = r] V h c l^h b + 0(l/m 2 ), 07^756 = r] A h c ~f^ 5 h b + 0(l/m 2 ). 



(3.28) 



They are given by a perturbative series, in which, as above, we only keep the BLM-type 
/3Ja™ +1 terms (z = m c /m b ): 



Tj V ,A = 1 + C F 



71 



i + £(-/? )X ,A <M 



n=l 



with m |a 



3 
2 



3 1 + 2 

4 1 - z 



In z, r A {z) 



-2 In 2; . 

4 1-2 



(3.29) 



(3.30) 



The resummation of (3q terms for the 77's was discussed in some detail in || |Tl|] and can 
equally easily be implemented within our dispersion technique. For completeness, we collect 
the necessary formulas in Appendix B. The results are summarized in Table H, where we 



give the leading-order coefficients d\ ,A and the resummed enhancement factors 



M% A [a, -f3 a s (m b )] = 1 + £ {-(3 ) n d n a n s {m b ) . 

n=l 



(3.31) 



To make an explicit comparison to inclusive decays possible, we have presented the results 
in the form of an expansion in a s (m b ) rather than in a s (y/m b m c ) which is more natural in 
exclusive decays. Because of this, our coefficients are related to the ones given in Ref. || 



by di' A (fi = y/m b m c ) = di' A (fi = m b ) + \wm c jm b . Since the product of a s and M ro is 
scale- independent, when a one-loop running coupling is used, we have 



a s (m b 



M^ A [z,-p Q a s { v fm b ~m~ c )] = , 

a s {y/m b m c ) 



M^ A [z,-p a s (m b )}. 



(3.32) 



Thus, for example, with a s (m b ) / a s ( ^/m b m c ) = 0.82 and M A [0.3, —/3oa s (m b )} = 1.96 from 
Table || we get M^[0.3, — j3oa s (^/m b m c )} = 1.59. Note that this number is rather large and 
indicates that the higher-order corrections are important in the axial channel. 
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m c /m b 


d\ 


d? 




M v 

oo 




M A 

oo 


0.1 


2.47 


1 


63 


2 


07 ±0 


05 


2.69 ± 1.09 


0.2 


1.78 


3 


32 


1 


50 ±0 


07 


2.37 ±0.42 


0.3 


1.36 


2 


56 


1 


30 ±0 


06 


1.96 ±0.16 


0.4 


1.08 


2 


17 


1 


20 ±0 


05 


1.76 ±0.09 


0.5 


0.86 


1 


90 


1 


14 ±0 


04 


1.64 ±0.06 


0.6 


0.68 


1 


70 


1 


09 ±0 


03 


1.56 ±0.05 


0.7 


0.52 


1 


53 


1 


06 ±0 


03 


1.50 ±0.04 


0.8 


0.39 


1 


39 


1 


04 ±0 


03 


1.45 ±0.03 


0.9 


0.27 


1 


27 


1 


02 ±0 


02 


1.41 ±0.03 


1 


0.17 


1 


17 




1.00 




1.37 ±0.03 



Table 3: The lowest-order coefficients d\' A and the resummed series M^ A for t]v,a in the 
MS scheme as functions of the ratio of pole masses m c /mb- Note that the expressions 
diverge for m c — > 0. For m c = rrib, t]v = 1 due to charge conservation. Input parameter: 
-/3 a s (m b ) = 0.14. 



Where a comparison is possible, our results agree with the values obtained in ||11|| . We 
obtain 

i] A = 0.943 ± 0.005 ± 0.010 ± 0.001, (3.33) 

where the first error gives the estimated renormalon uncertainty^, the second one the uncer- 
tainty coming from a s (mz), and the third one the uncertainty in the input quark masses. 
In the limit m c — > rrib the inclusive decay rate thus equals (in perturbation theory) 

r( B X c eP) = m^Z^l {i _ Cf ^ \ M*(a = 1) j , (3.34) 

which implies the relation 

M h -\a = l)^Mi{a = l) (3.35) 

and provides a non-trivial check of our calculation. Comparing the corresponding entries in 
Table [I] and Table |3] we indeed find agreement. 



4 Determination of 1146 



6 We note that this uncertainty, although of order Aq CD /m 2 , is numerically smaller than the estimate of 
explicit 1/m 2 corrections, see Sec. 4.2, which justifies the addition of these explicit corrections, disregarding 
their potential ambiguity connected with renormalons. 
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4.1 Theoretical input parameters 



The main parameters we need in order to determine \V c b\ are the b and the c quark masses. 
At present there seems to be no general consensus on their values, and the existing estimates 
are often controversial. The masses extracted from the spectroscopy of bb and cc mesons are 
usually given with very small errors, see e.g. Ref. ||25|| . However, the actual uncertainty is in 



this case hidden in their relation to the running masses at a certain hard scale, which we need 
in this paper. Thus, we prefer to rely on a less accurate (as far as numbers are concerned) 
direct determination of the MS b quark mass from QCD sum rules for lowest moments of 
e + e~ annihilation to heavy quarks [26, 27|. Due to the smaller scales involved, such estimates 
are less reliable for the c quark mass, so that we prefer to fix m c by a different method, see 
below. In this paper we use the following value for the MS b quark mass |28|, [29IPT: 



m b (m b ) = (4.23 ± 0.05) GeV. 



(4.1) 



We point out that m b was chosen as renormalization point mainly in order to conform to the 
standard choice in the literature, in the same way as a s is usually normalized by its value 
at the Z boson mass. Actually in our approach the question of the "natural" scale does not 
appear, at least to the extent that the approximation of summing higher order corrections 
due to vacuum polarization is good: all results are explicitly scale-invariant for a one-loop 
running coupling. In finite order perturbative calculations it may be more appropriate to 
choose a lower renormalization scale in order to minimize the size of uncalculated higher 
order terms. From ft4.1|) we get the pole mass0 



(5.05 ± 0.06) GeV, 



(4.2) 



where all radiative corrections of type /3 
quark pole mass were proposed in Refs. 
0- 



n n+ 



a 



1 are resummed. Very similar values for the b 



30. |2pl and are also indicated by lattice calculations 



In order to fix the c quark mass, we make use of the fact that the difference between 
the pole masses of two heavy quarks is free from many ambiguities intrinsic in the mass 
parameters themselves and can be determined to a good accuracy from the expansion 



m b — m c = ijib — mo + 



1 



Ai + 3A 2 +0(a s /m, 1/m 2 ), 



(4.3) 



7 The number given in ]28| , |2 
one-loop accuracy used in Refs 




literally corresponds to the Euclidian mass ml(p = —ml), but to the 
the difference between the Euclidian and the MS mass is negligible. 
8 Recall that any numerical value of the pole quark mass implies its proper definition. Our central value 
corresponds to the principal value prescription to sum the perturbative series that relates the pole to the 
short-distance MS mass, the error comes from the 50MeV uncertainty in mb(rh~b)- The freedom in choosing 
the summation prescription results in an additional uncertainty of mb of order 100 MeV, which is exactly 
cancelled by a corresponding uncertainty in the decay rate, see the discussion in Sec. 3.3. 
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where m# and mp are the B and D meson masses, respectively; A2 is given by 

A 2 ~ ~(m%. -m 2 B ) ~ 0.12 GeV 2 , (4.4) 

and — Ai/ (2m b ) is the kinetic energy of a heavy quark inside a B meson. For Ai an estimate 
is available from QCD sum rules [32|, Ai = —(0.6 ± 0.1) GeV 2 , which is however strongly 



correlated with the value of A = mg — m b and thus with the value of the b quark pole mass. 
The estimate quoted in |32| is the average of Ai = —0.5 GeV 2 and —0.7 GeV 2 obtained for 



A = 400 MeV and 500 MeV, respectively. The resummation of (3$ radiative corrections 
leads to a significant increase of the value of the b quark pole mass, and thus to a much 
lower value of A of order (200-300) MeV. Thus in principle the value for Ai to be used in our 
analysis should be obtained from a QCD sum rule using a small value of A and including 
the resummation of running coupling effects, which is not available. We have calculated the 
BLM-type a 2 s correction to the simpler sum rule for the leptonic decay constant fs in the 



static limit ||33|| , which also enters the sum rule for Ai as normalization factor, and found 
that this correction is very small (in other words: the BLM scale coincides with the "naive" 
hard scale). This may be accidental, however, and rather indicate that other corrections are 



important. We have also checked that the sum rule for Ai derived in |32[ becomes much 
less sensitive to the value of A if it is small, and even with arbitrary A it is not possible to 
push — Ai below (0.25-0.30) GeV 2 . Lacking a BLM-improved sum rule for Ai, we consider 
— (0.5 ± 0.2) GeV 2 as a fair estimate. With this value, one obtains^] from ( |4.3j ) 

m b - m c = (3.43 ± 0.04) GeV , (4.5) 

assuming that a s /m and 1/m 2 corrections in (|4.3| ) are negligible. Combining this result with 
the b quark pole mass in ( |4.2| ) we get the c quark pole mass 

m c = (1.62 ± 0.07) GeV (4.6) 

and the running mass 

m c (m c ) = (1.29 ± 0.06) GeV . (4.7) 



This value is consistent with determinations from QCD sum rules p6| , [27J . 

For completeness, we also give the corresponding values of the "one-loop pole masses" 
defined as 

mQ = m 6jC (m 6 , c )(l + C F °^^-\ : (4.8) 



7T 



m W = (1.50 ±0.06) GeV, 



m « = (4.63 ±0.05) GeV. (4.9) 



9 Using our technique it is possible to estimate the uncertainty of the relation ( |4.5| ) that is due to 
infra-red contributions suppressed by three powers of the quark mass (as mentioned above, the 1/m 2 renor- 
malon uncertainty is absent, at least to our approximation). A simple calculation yields <5(to{, — m c ) ~ 
m c /(47r|/3 |)exp[5/2- 3/(2|/3 K(m c ))], which is of order (5-8) MeV. 
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Our values are in agreement with those quoted in Ref. [pOfl . 



4.2 Extraction of \V c b\ 

Our discussion was so far restricted to the perturbative corrections to the b quark decay 
rate. In order to extract \V c b\ from the experimental data, we now have to specify the non- 
perturbative corrections, which are suppressed by two powers of the heavy quark masses. 
Summarizing the results of Refs. || |32], |34| we write 



T(B -> X c eu) = Y{B^ X c eu) pert 1 + — (4.10) 




with 5 NP = —(1.05 ± 0.10) GeV 2 , where the error comes from the uncertainty in Ai, cf. 
Sec. 4.1. 

As for the exclusive decays, the experimentally interesting quantity is the differential 
decay rate at zero recoil of the final state meson, which depends on the form factor 

:f(i) = tm(i + <W). (4.ii) 

The short-distance correction tja was already discussed in Sec. 3.5; numerical values are given 
there and in Table |6| below. The non-perturbative correction 8i/ m 2 was estimated in Refs. || 
as -(5.5 ± 2.5)% using X 1 = -0.4 GeV 2 .|3 

As experimental input we use the world average B° lifetime t#o = 1.5 ps |36|, the most 
recent measurement of the B° semileptonic branching ratio BgL — (10.9 ± 1.3)% [j37[] , and 



|KfejF(l)| = 0.0354 ± 0.0027 p8[ , where we have rescaled the latter value to be compatible 
with t b o = 1.5 ps. 

Our results are summarized in Fig. [|and in Tables [| and where we give \ V c f,\ as function 
of the running b quark mass. The determination with one-loop radiative corrections is plotted 
in Fig. f|(a) and using resummation in Fig. ^(b). The solid line represents the result from 
the MS, the long dashes from the OS calculation. The short dashes give \V c b\ from exclusive 
decays. The shaded areas illustrate the range of b quark mass values from Eq. (|4.1| ). The 
curves in Fig. |](b) are obtained for Ai = —0.5 GeV 2 , which corresponds roughly to a fixed 
difference of the pole masses of m& — m c = 3.43 GeV. 

The two sets of curves in Fig. f|(a) are obtained with two different ways to specify the 
c quark mass. The first way is to use exactly the same short-distance c quark mass as 
in Fig. f|(b): the curves labelled (i) and (ii) are obtained using the constraint ( f4.3[ ) with 
resummed pole masses, i.e. fixing m& — m c = 3.43 GeV, and calculating and m c (m c ), 
respectively, for a given value of m^mj). Thus, since all short-distance input parameters are 
chosen in precisely the same way as in Fig. §(b), the difference between the predictions for 
\V c b\ is entirely the effect of resummation. 



10 



An earlier estimate f| was -(8.5 ± 3)% for Ai = -0.54 GeV 2 . 
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m b (m b ) 
m b 

m c (m c ) m c 



4.0 
4.80 



-(9) 
-(9) 



4.1 
4.91 



-(9) 
-(9) 



4.2 
5.02 



+(9) 
-(9) 



4.3 
5.13 



+(9) 
-(9) 



4.4 
5.24 



+(9) 
-(9) 



1.18 
1.22 
1.26 
1.30 
1.34 



1 47~ (6) 

1 52~ (5) 
i - oz +(o) 

1 58~ (4) 
1 63~ (4) 
1.681$ 



-(3) 
+(3) 

-(3) 
-(3) 

-(3) 
+(3) 

0-048;$ 



0.045 



0.046 
0.047 



0.049 



-(3) 
+(3) 



0.043;$ 
0.044;$ 
0.045;$ 
0.046;$ 



0.039 



o.o4o; 

0.041 



-(2) 
-(2) 

'(2) 
'(2) 

-(2) 



o.o4i;$ 



0.042 



+ (2) 

(3) 
(2) 

(3) 
+ (2) 



0.037 



0.037; 

0.038 



0.039 



0.039 



'(2) 
"(2) 

'(2) 
"(2) 

-(2) 
"(2) 

"(2) 
+ (2) 

'(2) 
+ (2) 



0.034 



0.035 
0.035 



0.036 



0.037 



-(2) 

-(2) 

"(2) 
"(2) 

-(2) 
"(2) 

"(2) 
+ (2) 

"(2) 
+ (2) 



Table 4: jV^I calculated in the MS scheme for a wide range of quark masses. The central 
value is obtained for a s (m z ) = 0.117, the upper (lower) number in brackets gives the shift 
of the last digit if a s is put to 0.123 (0.111). The masses are given in GeV. Experimental 
input: B SL = 10.9%, t b o = 1.5 ps. 



Ax [GeV 2 ] 


m b (m b ) 
m b 

m b - m c 


4.0 

4.80-$ 


4.1 

4.91-$ 


4.2 

5.02-$ 


4.3 

r 1 o + (9) 
0.±0_( 9 ) 


4.4 

5.24-$ 


-0.3 


3.39 


0.044i;$ 


0.0426;$ 


0.0414;$ 


0.0402;$ 


0.0392;$ 


-0.4 


3.41 


0.0438;$ 


0.0423;$ 


o.o4ii;$ 


0.0399;$ 


0.0389;$ 


-0.5 


3.43 


0.0434;$ 


0.0420;$ 


0.0408;$ 


0.0397;$ 


0.0387;$ 


-0.6 


3.46 


0.043i;$ 


0.0417;$ 


0.0405;$ 


0.0394;$ 


0.0384;$ 


-0.7 


3.48 


0.0427;$ 


0.0414;$ 


0.0402;$ 


0.039i;$ 


0.0382;$ 



Table 5: \V cb \ calculated in the MS scheme in dependence on the b quark mass and Ai. 
Notations and experimental input as in the previous table. 
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Figure 4: The value of \V c b\ extracted from the inclusive B meson semileptonic decay 
rate to one-loop accuracy (a) [the curves (i) to (iv) are explained in the text] and after 
resumming (3q radiative corrections (b) as a function of the MS b quark mass for fixed 
Ai = —0.5 GeV 2 . The solid and long-dashed curves show the predictions obtained by using 
the MS and OS scheme, respectively. The central value coming from exclusive decays is shown 
by short dashes and the shaded area gives the interval of b quark mass values suggested by 



QCD sum rules, Eq. Q4.1|). Experimental input: t b o = 1.5 ps, B SL = 10.9%, a s (m z ) = 0.117. 



Although this choice is presumably the most clear way to show the effect of resummation 
of running coupling effects, it also may be slightly misleading. Indeed, the large difference in 
values of \V c b\ between the one- loop and resummed formulas is mainly due to the fact that 
the one-loop pole masses defined in this way do not satisfy Eq. (|4.3|): -m['' = 3.13 GeV. 
Thus we also try another choice, calculating the c quark mass by enforcing the constraint 
( [4.3|) expressed in terms of one-loop pole masses, so that — = 3.43 GeV. The 
results are shown in the form of curves (iii) and (iv). This choice is less instructive as far as 
the comparison between one-loop and resummed results is concerned, but is probably more 
attractive phenomenologically. On the other hand, we note that with this choice the value of 
the short- distance c quark mass becomes very low, m c (m c ) = 0.98 GeV [m^' = 1.18 GeV] for 
M&(m&) = 4.23 GeV, which is hardly consistent with the QCD sum rules for the charmonium 
system. The large difference in the resulting values for \V c b\ shows the dilemma all strict 
one-loop calculations are inevitably confronted with: it is impossible to relate the three 
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independently determined input parameters m^m^), m c (m c ) and Ai to each other within 
the errors by one-loop equations, although all of them provide valid phenomeno logical input. 
It is only after inclusion of higher order perturbative corrections in Eq. (|4.3| ) in form of the 
resummed pole masses that the three values appear to be consistent with each other, and 
it may be considered as a serious argument in favour of BLM-improved perturbation theory 
in semileptonic inclusive decays that the central value (but not the error bars, see below) of 
\V cb \ is independent of the choice of a particular subset. 

The effect of resummation is clearly visible in Fig. f|(b): first, we find a considerably 
reduced scheme dependence, i.e. the difference between the solid and the long-dashed curves 
is much smaller in Fig. f|(b) than in Fig. f|(a). Secondly, we observe good agreement between 
\V c b\ from exclusive and inclusive decays obtained with resummation, which otherwise is only 
achieved for either an unreasonably large b quark mass or a very small c quark mass. 

It has been proposed in Ref. that in inclusive decays the dependence on the quark 
masses is significantly reduced, if the charm and bottom masses are not varied independently, 
but related to each other by Eq. ([Of). To study this question, we give tables of numerical 
values for |V^|, choosing as independent input parameters either the running b and c quark 
masses (Table ^), or the running b quark mass and Ai, which specifies the difference between 
the pole masses (Table H). It is clearly seen that, although the central values for \V cb \ are 
nearly the same, the latter choice is preferable, since the inclusive decay rate is very sensitive 
to the mass difference — m c , and already a very modest accuracy in Ai in fact constrains 
nib—vn c more precisely[]]than any direct determination. We conclude that Ai (or equivalently 
mj — m c ) is a better theoretical input parameter than the c quark mass itself, in agreement 
with the discussion in Ref. [|35[| . In addition we find that the dependence on a s {m z ) is 



strongly reduced, too. We emphasize, however, that our choice of input parameters only 
serves to reduce the sensitivity of \V cb \ on uncertainties in the input parameters and that the 
central value is independent of that choice. 

It might also be useful to compare the results obtained from the resummed formulas with 
those obtained in the usual BLM approximation, where of the whole series of corrections 
generated by the running of the coupling only the a^Po term is taken into account. Again, 
the major subtlety comes from the necessity to specify the c quark mass. It turns out that 
the perturbative series that relates the c quark pole mass to the MS mass starts to diverge 



already in second order JTJ]. Because of the divergence, it is not justified to cut the series 
at that order, although in this case the BLM correction gives an excellent approximation 
to the exact two-loop result. For example, starting from m c (m c ) = 1.26 GeV one gets for 
the BLM pole mass = 1.75 GeV, which is significantly larger than the resummed 

result m c = 1.58 GeV. In other words, the BLM approximation underestimates the scale 
of the coupling and thus overestimates the radiative correction. A comparison between the 



11 Even if we abandoned the determination of Ai in Ref. f3l| completely and only put the constraint 
-0.7 GeV 2 < > 
at most 0.0014 



0.7 GeV 2 < Ai < OGeV 2 , which corresponds to 3.34 GeV < m b - m c < 3.48 GeV, \V cb \ would change by 
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BLM and the resummed results is however possible if one starts from the MS b quark mass, 
calculates the b quark pole mass in the BLM approximation, and then fixes the c quark pole 
mass from the heavy-quark expansion (}4.3|) . Taking for definiteness m^mj) = 4.23 GeV and 
Ai = —0.5 GeV 2 , we find in this way in the BLM approximation \V c b\ = 0.0381 (\V c b\ = 
0.0424) in the MS (OS) scheme, respectively, compared with |V cfe | = 0.0404 (\V cb \ = 0.0424) 
with resummed formulas. The perfect agreement between the BLM and resummed formulas 
in the OS scheme is probably accidental. In the MS (OS) scheme the BLM result is very 
close to the one-loop calculation, see Fig. [|. In both cases, the numerical effect of higher 
order corrections is small and of the order of 5%. We wish to emphasize once more, however, 
that this comparison is only possible if the c quark mass is obtained in a very particular way 
and it is only after resummation that we are able to get a selfconsistent description in terms 
of short-distance parameters. 

From the combined evidence of Fig. |] and Tables |] and [5] we extract the following results: 

\V cb \ m = (0.0404 ±0.0006 ±0.0006 ±0.0006) (7^7) ' ' , (4.12) 



VlO.9%/ V t b o 

\V cb \ os = (0.0424 ±0.0004 ±0.0007 ±0.0005) (7^7) ' (h!L&\ ' , (4.13) 



VlO.9% 



Tbo 



where the first error comes from the uncertainty in a s (mz) = 0.117 ± 0.006, the second one 
from the uncertainty of the b quark mass, Eq. (|4.1|) , and the third one from the uncertainty 
in Ai = —(0.5 ± 0.2) GeV 2 . For the reader's convenience, we collect all other parameters in 
Table § and Table 0. 

Comparing Eqs. ( f4.12| ) and ( |4. 13| ) , it is clear that the uncertainty of our calculation 



is dominated by scheme dependence, which is simply the effect of higher-order radiative 
corrections not related to the running of the QCD coupling and which we miss in our ap- 
proximation. A different way to estimate unknown higher-order corrections is to consider 
the scale dependence of the results in the MS scheme. This scale dependence is due to 
two-loop running effects in a s and terms of order (a s Ci?) 2 , which are beyond the accuracy of 
our approach. We have checked that after renormalization group improvement of ln(/i 2 /m 2 ) 
corrections (i.e. using the exponentiated form of the b quark mass scale dependence as in 
Eq. ( |3.20|) ) the variation of \V cb \ with the renormalization scale \x within m b /2 < \x < 2m b is 



of order ±0.001, i.e. of the same order as the difference between the OS and MS calculations. 
This may be taken as an estimate of the accuracy of the resummation, beyond which an 
explicit calculation of higher-order corrections is necessary. Combining the errors, we get 

(r B o/1.5ps) 1/2 \Vcb\ind = 0.041 ± 0.002 ± 0.002, (4.14) 

which is our final result. The first error shows the theoretical uncertainty and the second 
error comes from the experimental semileptonic branching ratio. The theoretical error is 
dominated by the uncalculated exact a 2 s correction to the decay rate. All large corrections 
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a s {m z ) 




m b [GeV] 


m c [GeV] 


a s (m b ) 


a s (m c ) 


M b ^ c [a, -/#V(m 6 )] 


VA 


0.111 


223 


4.96 


1.52 


0.190 


0.291 


2.07 


0.954 


0.117 


308 


5.05 


1.62 


0.208 


0.333 


2.07 


0.943 


0.123 


411 


5.15 


1.72 


0.228 


0.384 


2.02 


0.938 



Table 6: Input parameters in the calculation of the resummed B decay rate in the OS 
scheme and values of t]a entering the exclusive decay rate. 



a s (m z ) 




a s (m b ) 


a s (m c ) 


Ml(m b ) 


KM) 




(m 6 )] 


0.111 


0.287 


0.199 


0.326 


2.04 


1.82 


2.06 


2.18 


0.117 


0.306 


0.220 


0.383 


2.08 


1.55 


1.97 


2.40 


0.123 


0.329 


0.244 


0.451 


2.09 


1.22 


1.82 


2.68 



Table 7: Input parameters in the calculation of the resummed B decay rate in the MS 
scheme. 



coming from the running of the strong coupling either cancel after proper treatment of the 
infra-red regions or are cast into the redefinition of the scale in the coupling or, equivalently, 
into the "M-factors". The numerical significance of these M-factors depends on the way 



how the input parameters are chosen, and can be minimized by using the constraint ( |-1.3j ) 
following from the heavy quark expansion. 
In turn, we get from the exclusive decays 

(r B o/1.5ps) 1/2 \V cb \ exd = 0.040 ± 0.001 ± 0.003, (4.15) 

where the first error is the theoretical uncertainty^ and the second one the experimental 
error of the decay rate. 

At present the experimental errors are roughly the same for both the exclusive and the 
inclusive determinations. Actually the two approaches are complementary to each other: the 
inclusive decays offer a better opportunity to reduce the experimental errors, whereas the 
exclusive decay rates are inevitably plagued with large statistical errors from a measurement 
near the edge of phase space. On the other hand, in inclusive decays, the theoretical predic- 
tions are more sensitive to errors in the quark masses. It is encouraging that even now, with 
moderate experimental accuracy and an improvable accuracy of the theoretical input, both 
methods lead to very similar results. 



12 The theoretical error indicated in (4.15) does not include neither possible perturbative corrections, not 
related to running of the coupling, nor the uncertainty in Ai which contributes to the overall 1/m 2 correction. 
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ml • (n\ 

lhis paper 1 ' 


Ref. 39 W 


Ref. 40 (b) 


Ref. 18| (fe) 


Ref. [35] (a) 


Scheme 


MS OS 


OS 


OS 


MS 


OS 


m c [GeV] 


1.26^ 1.62^ 




1.57^ 


1.35^ 


1.3^ 


m b [GeV] 


4.23 (c) 5.05 (a!) 




4.96 (d) 


4.6 (c) 




\v cb \ 


0.040(1) 0.042(1) 


0.046(8) 


« 0.042 


0.036(3) 


0.042 



nib from QCD sum rules, m c from Eg . (|4.3| ). ^ Running mass normalized at m. 

( b ) m c from £>(D — > XeP), m& from Eg. (|4 . 3[) . ( d ) Pole mass. 



Table 8: |V^,| from inclusive decays obtained in previous analyses. If necessary, the quoted 
numbers are rescaled to be compatible with r B o = 1.5 ps, Bsl = 10.9%. The quark masses 
are the central values used in the papers. 

5 Conclusions 

We have carried out a detailed analysis of the radiative corrections to inclusive semileptonic 
B decays that originate from the running of the strong coupling. Independent of the actual 
accuracy of this approximation, our results clearly indicate that the series of radiative cor- 
rections in the OS scheme diverges, starting already in low orders. Thus determinations of 
\V cb \ from inclusive decays using the pole masses of b and c quarks as input parameters are 
plagued with a numerically large uncertainty. The accuracy cannot be improved by choosing 
a different scheme or a lower scale in the running coupling. 

However, the problem is spurious and entirely due to the use of a bad input parameter - 
the pole mass — which imports infra-red contributions at the level of 0{l/m b ) corrections. 
It can be shown |15], [TfJ that these corrections are absent in the inclusive decay rates, and 
therefore using pole masses in the tree-level decay rate induces large radiative corrections 
of infra-red origin simply in order to cancel infra-red effects hidden in the definition of the 
mass parameter. 

We demonstrate that the behaviour of the perturbative series is indeed drastically im- 
proved by using the MS mass instead. The calculation in the OS scheme can be saved, 
if the pole mass is defined by a certain non-perturbative prescription in its relation to the 
short-distance mass (or some physical quantity from which it is determined), and if the same 
prescription is used to sum the series of radiative corrections to the decay widths. This es- 
sentially implies a rearrangement of radiative corrections in two pieces, hiding part of them 
in the tree-level phase space, and in practice involves considerable cancellations already in 
low orders. 

We carry out a detailed analysis of inclusive decay rates with resummation of corrections 
induced by one-loop running of the coupling and find a good agreement between values of 
\V cb \ extracted from inclusive and exclusive decays. The comparison of our results with earlier 
calculations is presented in Table |8|. In general, we find agreement within the errors. It should 
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be emphasized, however, that predictions for inclusive decays depend rather strongly on the 
quark masses, and effects of the resummation can be masked by using different input values. 
We believe that one advantage of our approach lies in using well-defined mass parameters, 
which can (in principle) be extracted from experimental data with high accuracy. 

The accuracy of our predictions is limited by the unknown accuracy of the resummation 
of BLM-type radiative corrections. The incompleteness of this procedure is at least partially 
indicated by the scheme dependence of the result for |V^|, which is of order 5%. To reduce 
this remaining error it will be necessary to incorporate exact a 2 s corrections to the decay 
rates. One could try to reduce radiative corrections by expressing the decay widths in terms 
of masses renormalized at smaller scales, and adjust their values in order to reproduce the 
bulk of available data on heavy hadrons in the framework of one-loop calculations. This 
approach can be phenomenologically successful inasmuch as the structure of higher-order 
radiative corrections is similar in various processes, which is natural to assume, but difficult 
to control theoretically. 

On the other hand, with a restriction to the level of accuracy of order 5% for |V^|, we find 
that the theoretical calculation can be justified, and only a moderate precision is required 
for input parameters such as quark masses. In particular, it is sufficient to know the running 
b quark mass to an accuracy of order (50-100) MeV, and the difference in pole masses of b 
and c quarks to ~ 50 MeV, which corresponds to the uncertainty of the kinetic energy of the 
b quark inside the B meson of order 0.2 GeV 2 . We think that it is possible to achieve this 
kind of accuracy, and in fact even to improve it. 

To summarize, we conclude that inclusive decays of B mesons provide a valid source of 
information on |V c f,| with a present theoretical accuracy of order 5%, which presumably can 
be improved in the future. 
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A Radiative corrections to inclusive decays 

In this appendix we give an explicit expression for the one-loop radiative correction to the to- 
tal inclusive decay width. We use the on-shell renormalization scheme, where for finite gluon 
mass the wave-function renormalization constant is given by (in dimensional regularization 
in D dimensions): 

Z 2F = l + C F ^|-^ I + 7i ; -ln47r + ln4-4-3 2 /+(^ 2 -2) \ny 
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3(y 2 - 2y - 4) + 



2J1 - A/y 



(A.l) 



For a massive quark's self-energy £ we impose the renormalization condition S fi (m 2 ) = 0. 
Throughout the appendix, we use the notations y = X 2 /ml and a = ml /ml. 

The function do(a, A 2 ) can be written as the sum of contributions of real and virtual 
gluon emission, corresponding to different imaginary parts of the diagrams in Fig. [TJ: 



f 1 (a)g (a)d (a } \ 2 ) = - D mrt (a, y) + 6(m b - m c - X)D brems (a, y) 



(A.2) 



(recall that fx is the tree-level phase space factor defined in ( |1.2|) and go (a) the one-loop 
correction for zero gluon mass). The subscripts accompanying the -D's specify the diagrams 
in Fig. from which they are obtained, so that 



D = D I + D II + Dm + Djjj + D IV . 



(A.3) 



A.l The decay b — ► ceu 

For a massive quark in the final state, we did not succeed in obtaining an analytic expression. 
Below we give the renormalized contributions of the single diagrams, in Feynman gauge, 
expressed in terms of at most two-dimensional integrals: 



n f(\\ 1 3y 3{A + 2y-y 2 ) 1 + y/l -4/y ^ 3 2 1\ 1 

£>/ = A(aW -1 - — + . In v + l-y --)\ny\, (A.4) 

I 4 8^1 - Ay l-Jl-Ay V8 2) j 



= D// + fx [a] 



D 



brems 
II 



In a 



dx fx(x) w(a, x, y) 



y^y/a 

a 2 + x 2 — 6ax — y(a + x) 



Ax 2 (x — a) 2 



(A.5) 



(A.6) 



(l-v^) 2 1 

rffi + D\fp = J dx J dzw(l, a, x) 





In a 



-2xy H f-4a + 2y + 2a 2 + y 2 + 2a 3 + 2ax 

2a v 



O O 

+2ay + 2xy — 4a x — 4a y + 2ax — ay — 3xy — Ax y + 2axy 



+ w(l, a, x) (3x + 2y) In 



1 + Jl - 4a/(l + a- x) 



l-7l-4a/(l + a-x) 2 2 
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x(-2-y + 2a 2 - Ax 2 + 2ax + ay - 3xy) In 1 + 4 ^ + — Jl - — 

1 - y/l - A/y 2a V V 

1 + a/1 — Aa/y y\nv 

x (2 + 2x + y - 2a 2 - Ax 2 -ay- 3xy) In v , + - — - (-2 + 2a 

l-yjl- Aa/y 2a 

— 2x — y + 2a 2 — 2a 3 + 4a; 2 — Aax + 2ay — 2a 2 x — a 2 y + Aax 2 + 3xy + 3axy) 
1 

— — (1 — a — a 2 + a 3 — 3a; 2 + 2a; 3 + y + Aax — 3ax 2 — 2ay + a 2 y — 2xy 



( 1 1 + J1-4AAA 

+ 2xy 2 + Ax 2 y - 2axy) In A - lny + . In v , 

V ^/l - AA/y 1 - ^ - AA/y) 

(A.7) 

Dj^" + = /dx I dz w V> x >*) w W> x >y) ( 2 + a - x + y + 2 z + x 2 -Az 2 

J J x(a — x) v 

— 3ax + az — a;y — — 3yz) H f 1 — x + y + a; 2 — 3z 2 + 2z 3 — 2aa; 

a — x v 

— ay + 4az — xy — 2yz + a x — a z — az — 2xz + 4yz — 2y z + axy 

iu(l,x,2;)iy(a,x,j/) + (l-z + a;) (a-x-j/) + 2xj/ 

+ axz — 2a;y2;) In ■ 



Drr* = I ^f 1 (a/x)w(l,x,y) x2{1+x2 An 6X * {1 + x)} , (A.9) 



-u>(l, a;, z) w(a, x, y) + (1 — z + x) (a — x — y) + 2 xy 

' (A.8) 

(1 ~^ )2 )— 

r f, I n I ti ■?/)( 1 T ■?; I 

4(1 -a;) 2 

a 

where we have used 

w(x,y,z) = (x 2 + y 2 + z 2 -2xy -2xz -2yz) 1/2 , A = 1 - z + az - x(l - z)z. (A.10) 

The integrals have been further evaluated numerically. In doing so it is convenient to use 
the following expansion for large y: 

D mrt V =°° £ K(a) + (n(a) lny] ^ + 0(lny/y 4 ) (A.ll) 

n=l 2/ 
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0.99 


«2 
«3 


4321 
432 

667133 
14700 


5tT 2 

9 

87T 2 

3 


4.53 
19.1 


4.57 
19.2 


4.62 
19.4 


4.71 
19.6 


4.89 
19.9 


5.25 
20.6 


5.89 
22.1 


6.96 


8.64 


10.83 



Table A: Coefficient-functions in Eq. (|A.11|) . 



with the coefficient functions 
2 



/iCi = - - (l + 55a 2 -55a 3 - a 5 + 30a 2 lna + 30a 3 lna) , 

1 5a 2 
f x ( 2 = — (a 2 - l) (97 - 648 a + 1042a 2 - 648 a 3 + 97 a 4 ) + f 9 - 16a + 9a 2 ) lna, 

- — (1 - a) f 1349 - 8731 a - 3817 a 2 + 5598 a 3 - 3817 a 4 - 8731 a 5 + 1349 a 6 ) 
105 v ' 

+ 40 a 2 (1 + a) (3 - 4 a + 3 a 2 ) In a, 

^ (1 - a) (-41 + 284 a + 3834 a 2 + 284 a 3 - 41 a 4 ) 



/iC; 



Mi 



2a 



+ In a (-80 - 135 a -a 3 + 30 a lna) . 



(A.12) 



The coefficients k 2 , k 3 were calculated numerically and are given in Table ^ for several values 
of a. 

A. 2 The decay b — ► ueu 

In the limit a — > 0, the integrals in Eqs. (|A.6|) -( [^5 ) can be done analytically, yielding: 

D virt = ^ ( 14 V + 18 °y 2 - 2540 ?/ - 513 ) ~^(y 4 - l8 y 2 + + 6 



432 
1 



24J1 - 4/?/ 



4|/ 3 - lly 2 - 62y + 132) In 



1 + Jl- 4/y 



1 - a/1 - 4/y 



i(, 4 -18, 2 + 16, + 6){ln| + ln(l + ^)j{ln| + ln(l-yi~l)j 



32 



+ 



y 2 (y 2 + 2y- 12) C^l L 



6 



v 



u 



K'-v'-s 



+ ^ (l2y 3 + 45y 2 - 172y - 120) Iny 

i-y 



(A.13) 



brems 



(259^ + 43y 2 + 133y - 1863) - — - (36y J + 477y 2 - 172y - 120) 

I — 



432 

+ — In 2 y + I- (y 4 - 18y 2 + lQy + 6) ( 3 In 2 y - 2vr 2 - 12 arctan 2 



12 



72 



+ 24 arctan J- - 1 x arctan ( - J— - 1 ] + 12 < L 2 
\y \2-VVV J I 



+ L 2 



+ !^ !/ 2 (a 2 + 2y-12) 





/ « 4 \ 4 \ 4 
x < 2 arctan | W 1 +2 arctan J 1 + arctam / 1 — 7r 



z-y \ y 



3 — 2/ V 2/ 



+ — (1% 3 + 26y 2 - 234?/ + 72) J- - ll arctan J - - 1 - arctan 
36 V y V y 



2-y V V 



- 1 



+ 



36 JA/y - 1 



(19j/ 4 - 38y 3 - 299y 2 + 678?/ - 468) arctan 



i-y 4 

3-3/ V 3/ 



(A. 14) 



B Radiative corrections to exclusive decays 

The /JJa™ +1 corrections to exclusive decays at zero recoil were considered in Ref. [11|, using 



a somewhat different approach. To convert to our technique, it suffices to observe that for 

is 



ultraviolet convergent quantities the invariant distribution function w introduced in ITT 
related to the discontinuity of the one-loop radiative correction r with finite gluon mass A 
(in the normalization of Eq. ( |3.29| )): 



ro(A ) = 4C^J ds 7Tx 2 ~ s€,{T = s/fl ) - (Ri: 



o 
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Using /i 2 = m c mb and the explicit expressions^ given in JTT|, we find (with y = A 2 /m 2 and 
z = m c /m b ): 



r o(y) 



-3y + 2yz - Az 2 - 3yz 2 y (1 - zf (3y + Ayz - 2z 2 + 3yz 2 



+ 



v>z* 16z 4 
3y 2 — 5y 2 z — 2yz 2 + Qyz 3 + 6z 4 + Qz 5 



lay 



lz\l - z) 



Yaz 



3y 3 - 5y 3 z - 8y 2 z 2 + \<oy 2 z 3 + Ayz 4 + Ayz 5 - 32z 6 , 1 + a/ 1 - Az2 /v 

^=^= In — 



i - Vi - ** 2 /v 



lQyz 4 (l - z)Jl-Az 2 /y 



Ay - IQy 2 + 5y 3 + 32z - Ayz + 8y 2 z - 3y 3 z 1 + Jl-A/y 

la = , (^x3.2 J 



+ 



l&y(l-z)y/l-i/y 1-y/l-A/y 

9y + 2yz + 2Az 2 + 9yz 2 _ 9y 2 - 7y 2 z - Qyz 2 - Qyz 3 + I8z 4 + I8z b 
24^ 24(1 - z)z 4 

y {9y + 2yz - 6z 2 + 2yz 2 - 12z 3 + 2yz 3 - 6z 4 + 9yz^ 



In z 



A8z 4 



lay 



9y 3 - 7y 3 z - 2Ay 2 z 2 + 8y 2 z 3 + I2yz 4 + AAyz 5 - 96z 6 , 1 

^=^= In — 



1 -4^ 2 /y 



1- Jl-Az 2 /y 



A8yz 4 (l- z)Jl-Az 2 /y 



-AAy - 8y 2 + 7y 3 + 96^ - 12yz + 2Ay 2 z - 9y 3 z , 1 + V 1 ~ 4 A/ 

In , . (-d.o) 



4%(1 - zWl - 4/y 



13 Note that in Eq. (63) of Ref. Q (1 + t/2) must be replaced by (1 + t/z). 
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